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A sum rule analysis of the hadronic r-decay data can be used to determine the low-energy constants L\ (fj,) and 
Csf(fj,). These constants are QCD chiral-order parameters, which appear at order p 4 and p 6 , respectively, in the 
chiral perturbation theory expansion of the V — A correlator. At order p 4 we obtain L\ (M P ) — — (5.22±0.06)-10 -3 . 
Including in the analysis the order p 6 contributions, we get L r w (M p ) = -(4.06 ± 0.39) ■ 10~ 3 and C"J 7 (M P ) = 
(4.89 ±0.19) • 10~ 3 GeV" 2 . 
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1. Introduction 

The fact that the r is the only lepton mas- 
sive enough to decay into hadrons makes it an 
excellent tool to study QCD, both perturbative 
and non-perturbative. The main idea under these 
studies is the use of the analyticity properties 
of the different two-point correlation functions, 
which are relevant for the dinamical description 
of the r hadronic width. 

As it is well known, analyticity of two-point 
correlation functions allow us to relate different 
regions of the ^-complex plane. Roughly speak- 
ing, one can relate in this way regions where we 
are able to compute analytically, either with Chi- 
ral Perturbation Theory (xPT) or with the short- 
distance Operator Product Expansion (OPE), 
with regions of the complex plane where we are 
not able to compute (except in the lattice) but 
that are experimentally accessible. This connec- 
tion can be used either to predict observables that 
we are not able to calculate "directly" or, in the 
other way around, to extract the value of QCD 
parameters that are not fixed theoretically. 

An excellent example of this last strategy is 
the accurate determination of the QCD cou- 
pling a s (M T ) [1-5] from the inclusive r decay 
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width into hadrons, which becomes the most pre- 
cise determination of a s (Mz) after QCD run- 
ning. Other parameters of the Standard Model 
that have been extracted from r physics are the 
strange quark mass and the Cabibbo-Kobayashi- 
Maskawa quark- mixing | V us | [6-10]. 

Non-perturbative QCD quantities can also be 
obtained from r-decay data. The fact that the 
spectral function of the r decay can be separated 
experimentally in its vector and axial- vector con- 
tributions allows us to study their difference, that 
is specially interesting because it vanishes in per- 
turbative QCD (in the chiral limit) and therefore 
it is a purely non-perturbative quantity. 

The r-decay measurement of this V—A spectral 
function has been used to perform [11-13] phe- 
nomenological tests of the so-called Weinberg sum 
rules (WSRs) [14], to compute the electromag- 
netic mass difference between the charged and 
neutral pions [12], to determine several QCD vac- 
uum condensates [15, 16] and also to determine 
the AI = 3/2 contribution of the AS = 1 four- 
quark operators Qj and Qg to e' K / Sk, in the chi- 
ral limit [17]. 

Using xPT [18-20], the hadronic r-decay data 
can also be related to order parameters of the 
spontaneous chiral symmetry breaking (S^SB) of 
QCD. xPT is the effective field theory of QCD 
at very low energies that describes the physics 
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of the SxSB Nambu-Goldstone bosons through 
an expansion in external momenta and quark 
masses, with coefficients that are order param- 
eters of SxSB. At lowest order (LO), i.e. 0{p 2 ), 
all low-energy observables are described in terms 
of the pion decay constant f v ~ 92.4 MeV and 
the light quark condensate. At 0(p A ) : the SU(3) 
%PT Lagrangian contains 12 low-energy con- 
stants (LECs), Li=i,...,io and H\p [20], whereas 
at 0(p e ) we have 94 (23) additional parame- 
ters Cj=i ... 94 {CJl 1 ... 23) in the even (odd) in- 
trinsic parity sector [21]. These LECs are not 
fixed by symmetry requirements alone and have 
to be determined phenomenologically or using 
non-perturbative techniques. The Lj couplings 
have been determined in the past to an accept- 
able accuracy (a recent compilation can be found 
in ref. [22]), but much less well determined are 
the 0(p 6 ) couplings Cj. 

There has been a lot of recent activity to de- 
termine the chiral LECs from theory, using as 
much as possible QCD information [23-32]. This 
strong effort is motivated by the precision re- 
quired in present phenomenological applications, 
which makes necessary to include corrections of 
0(p 6 ) where the huge number of unknown cou- 
plings are the major source of theoretical uncer- 
tainty. 

Here we present an accurate determination of 
the ;\TT couplings L 10 and C$7 [33], using the 
most recent experimental data on hadronic r de- 
cays [34]. Previous work on L\o using r-decay 
data can be found in refs. [12,13,15,35]. Our 
analysis is the first one which includes the known 
two-loop xPT contributions and, therefore, pro- 
vides also the 0(p 6 ) coupling C$7- 

We will first introduce the sum rule relations 
that we will use, then we will show our results 
and finally we will compare them with other re- 
cent analytic results and hadronic r data deter- 
minations. 



2. Theoretical framework: the sum rules 

The basic objects of the theoretical analysis are 
two-point correlation functions of the non-strange 
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Figure 1. Analytic structure of II(s). 

vector and axial-vector quark currents 
U^jiq) = i |d 4 2 ;e^(0|r(^(a ; )^(0)t)|0) 

= (^V+gV)ngJ 7 (g 2 ) + gV<i r (g 2 ) ) (1) 

where Jh denotes V^ d = u'y ll d and A^ d = ^7^75 d. 
In particular, we are interested in the difference 

II(s) = Tlf^y — ^ud^A ' an< ^ we w * n wor k m the 
isospin limit (m u = rrid) where n^y(g 2 ) = 0. 
The analytic behaviour of this correlator is shown 
in Fig[IJ together with the complex circuit that 
we will use to apply the Cauchy's theorem. 

As we are interested in relating the xPT do- 
main (very low energies) with the r data, we mul- 
tiply this correlator by a weight function of the 
form l/s" (with n > 0). In this way we gener- 
ate a residue at s — when we integrate over 
the circuit of Fig. Q] and apply Cauchy's theo- 
rem. Taking into account the OPE associated 
with our correlator at large momenta and work- 
ing in particular with the cases n — 1,2, one gets 
the following sum rules (see ref. [33] for a careful 
derivation) : 

r°° ris 1 2 f 2 

-82$= / — imn( s ) = 4 + n(o), (2) 

Js th s 7T mi 
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160J?s £^I to n W = ^ + f(o) < 3) 

where the integrations start at the threshold 
Sth = 4m 2 . These two relations represent the 
starting point of our work and define the effective 
parameters L°q and Cgy . Their interest stems 
from the fact that the l.h.s. can be extracted from 
the data (see Section [3]), while the r.h.s. can be 
rigourously calculated within xPT in terms of the 
LECs that we want to determine (see Section UJ. 

3. The r-data side 

We will use the recent ALEPH data on 
hadronic r decays [34] , that provide the most pre- 
cise measurement of the V — A spectral function. 
In the integrals of equations (J2J) and ((3|) we are 
forced to cut the integration at a finite value so, 
neglecting in this way the rest of the integral from 
so to infinity. The superconvergence properties of 
n(s) at large momenta imply a tiny contribution 
from the neglected range of integration, provided 
So is large enough. Nevertheless, this generates a 
theoretical error called quark-hadron duality vi- 
olation (Dvjf|. From the so-sensitivity of the ef- 
fective parameters one can assess the size of this 
error. 

In Fig. [21 we plot the value of L^o obtained for 
different values of so, with the one-sigma experi- 
mental error band, and we can see a quite stable 
result at so > 2 GeV 2 (solid lines). The weight 
function 1/s decreases the impact of the high- 
energy region, minimising the DV; the resulting 
integral appears then to be much better behaved 
than the sum rules with s n (n > 0) weights. 

There are some possible strategies to estimate 
the value of and his error. One is to give 
the predictions fixing so at the so-called "dual- 
ity points" , two points where the first and second 
WSRs [14] happen to be satisfied. In this way 
we get = -(6.50 ± 0.13) • 10" 3 , where the 
uncertainty covers the values obtained at the two 
"duality points" . If we assume that the integral 
(jSJ oscillates around his asymptotic value with 

3 Equivalently, we are assuming that the OPE is a good 
approximation for at any |s|=so, what is not expected 
to happen near the real axis, and that produces the DV. 




Figure 2. L|o(sq) and C|f (so) from different sum 
rules. For clarity, we do not include the error 
bands associated with the modified weights. 



decreasing oscillations and we perform an average 
between the maxima and minima of the oscilla- 
tions we get LIq = -(6.5 ± 0.2) • 10~ 3 . Another 
way of estimating the DV uses appropriate os- 
cillating functions defined in [36] which mimic 
the real quark-hadron oscillations above the data. 
These functions are defined such that they match 
the data at ~3 GeV 2 , go to zero with decreasing 
oscillations and satisfy the two WSRs. We find in 
this way L°q = -(6.50 ± 0.12) ■ 10" 3 , where the 
error spans the range generated by the different 
functions used. Finally we can take advantage of 
the WSRs to construct modified sum rules with 
weight factors w(s) proportional to (1 — s/sq), 
in order to suppress numerically the role of the 
suspect region around s ~ so [2]. Fig. [2] shows 
the results obtained with u>i(s) = (1 — s/sq)/s 
(dashed line) and u> 2 (s) = (1 — s/so) /s (dot- 
dashed line). These weights give rise to very sta- 
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ble results over a quite wide range of Sq values. 
One gets Lff = -(6.51 ±0.06) • 10" 3 using Wi(s) 
and iff = - (6.45 ± 0.06) • 10~ 3 using w 2 (s). Tak- 
ing into account all the previous discussion, we 
quote as our final result: 

L c $ = -(6.48 ±0.06) • 10~ 3 . (4) 

We have made a completely analogous analysis 
to determine Cf 7 . The results are shown in Fig. [2] 
The solid lines, obtained from Eq. ([3|), are much 
more stable than the corresponding results for 
L°q, due to the 1/s 2 factor in the integrand. The 
dashed and dot-dashed lines have been obtained 
with the modified weights W3(s) = ^ fl — and 

w 4 (s) =jt (l-^Y (l + 2^). The agreement 
among the different estimates is quite remarkable, 
and our final result is 

C$ = (8.18 ±0.14)- 10" 3 GcV" 2 . (5) 

4. The xPT side 

Using the results of ref. [37] to calculate within 
%PT the l.h.s. of equations (J2j) and ((3]) we get 

-8Lf a = -SL^+GtM 

+ G 6 0L (n) + G 6 1L (») + G e 2 L (ri, (6) 
16C 8 c 7 ff = H\ L 

+ 16C< 87 ( f i) + H 6 1L (ri + Hl L ( f i),(7) 

where the functions H™ L (ji) are correc- 

tions of order p m generated at the n-loop level, 
which explicit analytic form [33] is omitted for 
simplicity. 

Working at 0(p A ) the determination is 
straightforward and one gets 

L r w (fi = M p ) = -(5.22 ± 0.06) • 10~ 3 . (8) 

At order p 6 , the numerical relation is more in- 
volved because the small corrections G® L 1L (fJ>) 
contain some LECs that represent the main 
source of uncertainty for L\ . It is useful to clas- 
sify the 0(p e ) contributions through their order- 
ing within the 1/Nq expansion. The tree-level 
term Gq L (/x) contains the only 0(p e ) correction 
in the large-Afc limit, 4m 2 (Cg 1 -C I 12 -C S0 ); this 
correction is numerically small because of the m 2 



suppression and can be estimated with a mod- 
erate accuracy [27,28,37-39]. At NLO G% L (p) 
contributes with a term of the form rn 2 K (C§ 2 — 
C[ 3 — Cgi). In the absence of information about 
these LECs we will adopt the conservative range 
\Cl2-Cls-CU < \C r 61 -C r 12 -C r 80 \/3, which gener- 
ates the uncertainty that will dominate our final 
error on L\ . Also at this order in 1/Nc there 
is the one-loop correction G\ L {n), that is pro- 
portional to Lg which is better known [40]. Cal- 
culating the 1/Nq suppressed two- loop function 
G\ L (/i) and taking all these contributions into ac- 
count we finally get the wanted 0(p 6 ) result: 

L r w {M p ) = -(4.06 ± 0.04 i? « ± 0.39 LECs ) • l(T a 

= -(4.06 ±0.39) • 10~ 3 . (9) 

where the error has been split into its two main 
components. Repeating the same process with 
C 87 (where the only LEC involved is Lg) we get 

C 8 r 7 (M p ) = (4.89 ± 0.19) • 10~ 3 GeV~ 2 . (10) 

5. Summary 

Through a sum rule analysis, that only uses 
general properties of QCD and the measured V— A 
spectral function [34], we have determined the 
chiral LECs L\ Q (M p ) and C 87 (M p ) rather accu- 
rately, with a careful analysis of the theoretical 
uncertainties. 

There are other determinations of L10 from r 
data in the literature. Our result for agrees 
with [12,15,35], but our estimation includes a 
more careful assessment of the theoretical er- 
rors. The 3.2 a discrepancy between the esti- 
mation of ref. [13] and ours is caused by an un- 
derestimation of the systematic error associated 
with the duality-point approach used in that ref- 
erence. In ref. [35] C| 7 is also determined, in 
good agreement with our result. The extraction 
of L iq [p) from L^o nas only been done previously 
in ref. [12], at 0(p 4 ). 

Our determinations of L\ Q (M p ) and C 87 (Af p ) 
agree within errors with the large-iV^ estimates 
based on lowest-meson dominance [24,29,37,42], 
L w a -3/ 2 /(8Af 2 ) « -5.4 • 10" 3 and C 87 w 
7/ 2 /(32M£) « 5.3 • 10" 3 GeV~ 2 , and with the 
result of ref. [31] for Cg7, based on Pade approxi- 
mants. These predictions, however, are unable to 
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fix the scale dependence which is of higher-order 
in 1/Nc- More recently, the resonance chiral the- 
ory Lagrangian [29,43] has been used to analyse 
the correlator II(s) at NLO order in the 1/Nc ex- 
pansion. Matching the effective field theory de- 
scription with the short-distance QCD behaviour, 
the two LECs are determined, keeping full con- 
trol of their fi dependence. The theoretically pre- 
dicted values L r w {M p ) = -(4.4 ± 0.9) • 1CT 3 and 
Cl 7 (M p ) = (3.6±1.3)-lCT 3 GeV- 2 [32] are in per- 
fect agreement with our determinations, although 
less precise. A recent lattice estimate [44] finds 
L r w (M p ) = -(5.2±0.5)-lCT 3 at order p 4 , in good 
agreement with our result (HJ). 

Using the results of ref. [41], the SU(2) \PT 
LEC 1$ can be extracted from L\ (n). We find 
1 5 = 13.30 ± 0.11 at 0{p A ) and l 5 = 12.24 ± 0.21 
at 0(p 6 ). 

Recent analyses of the decay ir + — > l + vj at 
0(p 6 ) have provided accurate values for the com- 
binations Lg + Lio [39] and — [45], that can 
be combined with our results to get Lg(M p ) = 
(5.5 ± 0.4) • 10~ 3 and l e = 15.22 ± 0.39 to order 
p 6 , in perfect agreement with refs. [40,46]. 
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